Higher twist corrections to the sum rule for semileptonic B decay 



Changhao Jin 

School of Physics, University of Melbourne 
Victoria 3010, Australia 



Abstract 

The sum rule for charmless inclusive semileptonic i?-meson decays allows a 
theoretically clean and experimentally efficient determination of \Vub\- The 
leading twist contribution to the sum rule is known in QCD. We compute 
higher twist corrections to the sum rule using the heavy quark effective theory. 
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A new method to determine the Cabibbo-Kobayashi-Maskawa (CKM) matrix element 
\Vub\ has been proposed |l| that takes advantage of the sum rule for charmless inclusive 
semileptonic 5- meson decays B X^ii^e {i = e or fi). The sum rule establishes a clean 
relationship between |T4fe| and the observable 

/•I 1 dr - 

S= / d^uj^-^iB^XJu,) (1) 

with the kinematic variable = {q^ + \q\)/MB in the 5-meson rest frame, where q is the 
momentum transfer to the lepton pair and Mb denotes the B meson mass. Moreover, this 
method of determining \Vub\ has experimental virtue too. The kinematic variable ^„ is the 
most efficient discriminator between B — > X^^vi signal and B —>■ X^ii^i background. A 
majority of i? X^iui events have a value of C,u beyond the limit allowed for B X^ii^i 
decays with charm in the final state, > 1 — Mb/Mb = 0.65 with Mb being the D meson 
mass. Therefore, only a small extrapolation is needed to obtain S. 

The charmless inclusive semileptonic decay of the B meson is a light-cone dominated pro- 
cess. The light-cone expansion allows a rigorous and systematic ordering of nonperturbative 
QCD effects, providing an effective technique for a separation and classification of higher 
twist (HT) effects [0]. The leading term in this expansion gives the leading twist contribu- 
tion. HT contributions are contained in the light-cone expansion beyond the leading order. 
The sum rule at the leading twist order measures the bottomness carried by a -B meson. 
There are no perturbative QCD corrections to the sum rule. Thus the primary hadronic 
uncertainty and the potential uncertainty of perturbative QCD are eliminated, dramatically 
reducing the theoretical error on \Vub\- This inclusive method is to be contrasted with the 
determination of \Vub\ from the charmless inclusive semileptonic branching fraction of B 
mesons where the calculation of the total semileptonic decay rate is model dependent or 
assumes quark-hadron duality, there are uncertainties due to perturbative QCD corrections 
and, in addition, a larger extrapolation is necessary to extract the total rate if the kinematic 
cut on a certain observable, such as the charged-lepton energy or the invariant mass of the 
lepton pair, is applied for the suppression of 6 ^ c background. 

Only uncertainties due to HT effects remain in the sum rule. Including the HT contri- 
bution AHT, the sum rule reads 

Although HT contributions are expected to be suppressed by powers of Aqcj-,/M^ (Aqcd 
being the QCD scale), a quantitative estimate of them is indispensable for a complete 
understanding of remaining theoretical uncertainties in this determination of \Vub\- In this 
paper, we investigate HT effects on the sum rule for charmless inclusive semileptonic B 
decays using the heavy quark effective theory (HQET) [§-§1. 

Charmless inclusive semileptonic decays of the B meson are induced by the weak inter- 
actions. The differential decay rate to lowest order in the weak interactions is 

G% \V,ih\ , ,, (Pko d?k,, , ^ 
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Here E (P), Ei (ke), and E^, {k^) denote the energies (four-momentums) of the B meson, 
the charged lepton, and the antineutrino, respectively. The leptonic tensor for the lepton 
pair is completely determined by the standard electroweak theory since leptons do not have 
strong interactions: 

L^" = 2{e,k: + k^^k'; - g'^'ke ■ k, + le^''' afsk^k^^). (4) 

The hadronic tensor incorporates all nonperturbative QCD physics for the inclusive semilep- 
tonic B decay. It is summed over all hadronic final states and can be expressed in terms of 
a current commutator taken between the B meson states: 



~JdSe^'^y{B\[Uy),jl{0)]\B), (5) 



where jfj.{y) = u{ii)'^^{l — 75)&(?/) is the charged weak current for the h ^ u transi- 
tion. We adopt a covariant normalization for one-particle states, i.e., {B{P)\B{P')) = 
(27r)32p05(3)(p _p'). 

The most general hadronic tensor form that can be constructed is a linear combina- 
tion of P^Pu, P^iQu, (l^Pv^ (l^J.(lu, ^nvapP'^Q.^ ^ and g^y, with coefficients being scalar functions 
H^a(?7, q^) of the two independent Lorentz invariants, v = q ■ P/Mb and g^. However, the 
combination P^qv — q^Py does not contribute since L^^ {P^q^ — q^Py) = 0. Thus the hadronic 
tensor must take the form 

Equation (||) shows that VF^^ = W^^j,, so Wa, a = 1, . . . , 5 are real. The interesting physics 
describing the hadron structure and the strong interactions is wrapped up in the five dimen- 
sionless real structure functions Wa{i^, g^), a = 1, . . . , 5 for the unpolarized processes. 

In the following we will neglect the masses of the charged lepton and the w-quark. From 
Eqs. (H) and (^, we obtain the double differential decay rate for B Xuii'e in the rest 
frame of the B meson 

J2p ^2 IT/ |2 I 12 

1^1 iW,3q' + W,\ci\'), (7) 



where 



H = - (8) 



By integrating Eq. over g^, one gets the decay distribution of the kinematic variable 

dV r^il& , , d^Y 



d^u Jo d^udq^ 
Computing the current commutators one obtains from Eq. (j^) 

= -^{S^^^p - ze^^^p) J d'ye"i-y [9" A„(y)] {B \b{0)j^U{0,y)b{y)\B) , (10) 
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where S^^up = 9,iaguf3 + 9tif39ua - 9fi,ygaf3- In the above we have used 

{u{x), u{y)} = i{-f ■ d)iAu{x - y)U{x, y) (11) 

with the Wilson hnk 

U{x,y) = Vexp[tg, f dzf^A^^z)], (12) 

^u{y) = / d'ke-^'^yeik')S{k'), (13) 

where is the background gluon field and e{x) satisfies £(|a;|) = 1 and £:(— |x|) = —1. 

The matrix element {B\b{0)'y^U{0, y)h{y)\B) is the basic building block of the description 
of inclusive B decays in QCD. In general one can decompose it in the following form: 

(S|6(0)7^f/(0, y)h{y)\B) = 2[P^F{y\ y ■ P) + yPG{y\ y ■ P)l (14) 

where -F(y^, y ■ P) and G(2/^, y ■ P) are functions of the two independent Lorentz scalars, y"^ 
and y ■ P. The dominant part of the integrand in the hadronic tensor (^UJ) stems from the 
space-time region near the light cone, with the deviation from the light cone being of the 
order of the inverse large momentum ~ ~ 1/ -^00. The light-cone expansion 
of the functions -F(?/^, y ■ P) and Giy'^, y ■ P) m. powers of y"^ leads to 



{B\m^^U{Q,y)h{y)\B) 



n=0 n=0 



= 2 {P^ [T^'\y . P) + y'j^^'\y . P) + ...]+ / [g(^){y . P) + y^g«^){y ■ P) + ■■■]} . (15) 

The coefficients J^(2n+2)(^ . ^.^^ gC2^+^)(^y . p) in the light -cone expansion can be classified 
by twist. Following the notion of twist introduced by Jaffe and Ji [||, jF(^"+^)(y ■ P) has 
twist 2n + 2 and g^'^"'~^'^\y ■ P) has twist 2n-|-4, as from dimension analysis we know that the 
contribution of the former is suppressed by (Aqcd/^b)^", and the contribution of the latter 
is suppressed by (Aqcd/^b)^""''^- We will further discuss the non-local light-cone expansion 
of matrix elements below. 

The twist decomposition for the decay rate thus takes the form 

oo 

(ir = ^rfr(2«+2), (16) 

n=0 

where 

^r(2n+2) _ iKfel rf^'^U/('in+2) d'^^e /-,7\ 

" (27r)5E 2E, 2E, 

is the twist-2n + 2 contribution to the decay rate with 

TT/(2n+2) _ ^ TT/(2n+2) -^V^TT/{2n+2) • P°'q^ ^j^{2n+2) 

H B 
<hJly.rA2n+2) PfiQu + q,iPu T;f^(2n+2) 

' — Ml — ^^^^ 



X 



P^(y2)"^(2"+2)(^ • P) + /(2/')"-'e?('"+'^(2/ ■ P)l . (19) 
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The leading twist contribution to the sum rule @ results from JF^^) (y ■ P) of twist 2 and 
is known in QCD to be |T|] 

'd^u^^iB^XJue) = (20) 



which is a consequence of the conservation of the 6-quark vector current by the strong 
interactions. The next-to-leading twist contribution to the sum rule arises from JF^^) (y ■ P) 
and g^^\y ■ P) of twist 4. It can be obtained by integrating Eq. (0) over ^„ with the two 
relevant structure functions Wi^\v^q^) and W2'\h'^q^) of twist 4. 

We use the operator product expansion and the heavy quark effective theory to compute 
the twist-4 structure functions. The Wilson hnk is a gauge dependent operator. It is 
convenient to use the Fock-Schwinger gauge such that U (0, y) is unity. Since the h quark 
inside the B meson behaves as almost free due to its large mass, relative to which its binding 
to the light constituents is weak, one can extract the large space-time dependence 

6(^)=e--^-^6,(2/), (21) 

where mt is the 6-quark mass and v = P/Mb is the four- velocity of the B meson. This 
factorization makes clear why the large scale in matrix elements does not affect the relative 
size of terms in the light-cone expansion (plsl). The large scale hidden in matrix elements of 
6-quark operators is contained in an overall factor e-^^b^-y ^ gQ reduced matrix elements of the 
operators containing the rescaled operator involve only momenta of order Aqcd, which 
determine the relative size of terms in the light-cone expansion ([l5|), i.e., schematically 



oo /a2 \n 

{B\b{oh^Ky)\B) = e-'"^>'^-y{B\K{oh^K{y)\B) - e~'^^^-y x: (^^j • (22) 

The rescaled operator for a free 6-quark no longer depends on the space-time, so h{y) = 
e-'^^o^-ybiO). In this case all the coefficients jr(2"+2) . p) and ^(2n+4)(^.p) ^^e light-cone 
expansion (|15D vanish except that !F^'^\y ■ P) = e"*™''"'^, because the conservation of the 
6-quark vector current implies that {B\b{0)'y^b{0)\B) = 2P^. The leading-twist sum rule 
(pop is consistently reproduced in the free quark decay b —>■ uiu^. The conserved vector 
current 67^6 is not renormalized by the strong interactions. This explains why there are no 
perturbative QCD corrections to the sum rule. 

A Taylor expansion of the field in a gauge-covariant form relates the bilocal and local 
operators. This leads to an operator product expansion 



b{0)j%{y) = e-''^'^--yK{0)j\{y) = e"*'"''^-^ " " " vM^h'^k^'' " " " ^'^"^^.(O), 

(23) 



„=o 



where /c^ = iD^ = i{d^ — iggA^) and the symbol {■ ■ ■} means symmetrization with respect 
to the enclosed indices. Because of the weak dependence of the rescaled operator by{y) on y, 
we attempt to estimate the matrix element of the bilocal operator sandwiched between the 
B meson states with the truncated ^/-expansion in Eq. p3|). To obtain a twist-4 accuracy it 
suffices to keep only the first three terms 
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{Bm^%{y)\B) 



{B\KiO)^\{0)\B) 



+ i-t)y^{B\UO)^^tD^b,{0)\B) 



(24) 



In the heavy quark effective theory the QCD 6-quark field b{y) is related to its HQET 
counterpart h{y) by means of an expansion in powers of l/mb. 



Ky) 



-im,i,v-y 



1 + 



7p 

2mh 



+ 



'A 



QCD 



mt 



h{y). 



(25) 



The effective Lagrangian takes the form 



£hqet = hiv ■ Dh + h^^^h + h^-^^^^h + O 



2mh 



Arrih 



2 ; ' 



(26) 



where QsG^'^ = i[D^,D'^] is the gluon field-strength tensor. At the level of accuracy of the 
present discussion we take into account only the leading, l/rrib correction to the heavy quark 
limit rrib —>■ oo. Relating the matrix elements of the local operators in full QCD in Eq. (|24D 
to those in HQET, it follows that 



{B\b{0)^%{y)\B) = 2e~''^^''H 



1 - 2/ ■ ^^T^TT^^^ - iy ■ py 



3M, 



B 



+y^i-mhMBEh 



(27) 



where Eb = Kb + Gb and Kb and Gb are the dimensionless HQET parameters of order 
(AqcD/w^fc)^, which are often referred to by the alternate names Ai = —2mlKb and A2 = 
—2mlGb/3, defined as 



Comparing Eq. 



Ai 
A2 



2Mb 
1 



{B\h{iDyh\B), 
{B\hg,G^,a^'^h\B). 



12Mb 
with Eq. (^) yields 

J^^'\yP) = ^mlKbe-''^^^y, 
g^^\yP)=t'^mbMBEbe-'"'''^-y. 

We observe that the coefficients J-'^^\y ■ P) and Q^^\y ■ P) of the light-cone expansion (|I5|) 
are indeed of order Aq^^j-, as expected. Substituting Eqs. (|30D and (|3T|) in Eq. (|19[) and 
integrating by parts, we arrive at 



(28) 
(29) 

(30) 
(31) 
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+Ebe{q^ — mbV^)5{(f — 2mfcZ/ + ml) 
+mb{mb — v)Kbe[(^ — mbV^)5'{q^ — 2mfeZ/ + m^) 
+ — 2m;,z/ + ml)Ebe{q^ — mbV^)5\q^ — 2m6Z/ + m 
■1 
.2 

■1 



+ - mbV^)5'{q^ - 2m;,z/ + m^) 



+q^q,2MlEbe{q'' - mbv'')5'{q^ - 2m,z/ + 



-{Pt^Qu + qt,P^)mbMB [ - -KbX - Kbe{q^ - mbV^)5'{q^ - 2mbV + m; 



-2Ebe{q^ - mbv'^)5'{q^ - 2mbU + m' 



(32) 



where 6'(x) 



■£:S{x) and 



X 



dqf'dq^ 

Comparing Eq. (p^ ) with Eq. (|18|), we find 



— mbV^)5{q^ — 2mbZ/ + m^) 



(33) 



W['\u,q') 



16 11 

—mbMBi -mb{mb - y)KbX + Ebe{q^ - mbV^)6{q^ - 2m6Z/ + 



+ [mb{mb - u)Kb + (g^ - 2mfeZ/ + m^)^;,] 



Wt\^,q') 



xe{q^ — mbV^)5'{q^ — 2mbU + ml) 



16 



-KbX + 2{Kb + Eb)e{q'^ - mbV^)6\q^ - 2m,i/ + ml) 



16 

-—mlMlKb 



-X + £(g° - mbV°)6'{q^ - 2mfei/ + mf) 



32 



ym;,M|E5£(g° - mfet;°)(5'(g2 - 2m;,z/ + m^), 



16 9 , ^9 



-i^bX + {Kb + 2Eb)£(g° - mfei;°)y(g2 _ 2mbiy + m^ 



The twist-4 contribution to the sum rule can be obtained from Eqs. (H), (^, 
fi. The result is 



(34) 

(35) 
(36) 
(37) 

(38) 
and 



1 1 r/r(^) 



X 



IK 



1927r3 
80 m? 



304 



45 



76 

45' 



jr, + -Eb + 2-^Kb 



68 mf 



9 M| 
32 

" 15 M| 



26 mt 28 mt 112 mf 



3 M| 



3 M| 



15 M| 



(39) 
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This can serve as an estimate of HT contributions to the sum rule (0). 

For the numerical analysis, we need to know the values for the parameters involved. The 
HQET parameter A2 can be extracted from the B* — B mass splitting: A2 = (M|*-M|)/4 ~ 
0.12 GeV^, while Ai and mf, are less determined. For the purpose of estimation, we take 
Ai = -0.5 GeV^ and = 4.9 GeV 0. From Eq. (H), the HT correction to the sum 
rule (H) is then estimated to be AHT = 0.012. This quantitative study shows that HT 
corrections are at the expected level of ~ Aqq-[)/^b- 

In summary, we have elaborated a quantitative way of estimating HT contributions in 
inclusive B decays, which is based on the heavy quark effective theory. As an application we 
have calculated the twist-4 correction to the sum rule (0) for charmless inclusive semileptonic 
B decays. Using the sum rule, \Vub\ can be determined from a measurement of the weighted 
integral S. The error on \Vub\ due to HT corrections to the sum rule is estimated to be 
1%. Combining theoretical cleanliness and experimental efficiency, together with the better 
understanding of remaining theoretical uncertainties, the sum rule holds high promise of a 
precise and model- independent determination of \Vy,h\- 
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